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We consider the process of spin relaxation in the oscillating cantilever- 
driven adiabatic reversals technique in magnetic resonance force microscopy. 
We simulated the spin relaxation caused by thermal excitations of the high 
frequency cantilever modes in the region of the Rabi frequency of the spin sub- 



, system. The minimum relaxation time obtained in our simulations is greater 

m 

but of the same order of magnitude as one measured in recent experiments. 
We demonstrated that using a cantilever with nonuniform cross-sectional area 
Ch ' may significantly increase spin relaxation time. 

a ' 
cr. 

I. INTRODUCTION 

The method of oscillating cantilever-driven adiabatic reversals (OSCAR) is a promising 
method for spin detection in magnetic resonance force microscopy (MRFM) [1,2]. Just like 
the ordinary MRFM technique [3,4], the magnetic moment of a sample driven by a reso- 
nant microwave magnetic field interacts with a ferromagnetic particle. If the ferromagnetic 
particle is attached to the cantilever tip, the magnetic moment of the sample changes the 
cantilever vibrations. This change can be detected through measurement of the cantilever 
vibration parameters. In the OSCAR technique the cantilever is driven by an external force 
with a feedback loop designed to keep the amplitude constant. The cantilever vibrations 
together with the resonant microwave magnetic field cause a cyclic adiabatic inversion of the 
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magnetic moment of the sample. In turn, this cyclic inversion causes a shift of the cantilever 
vibration frequency, which is measured. 

The main challenge in the OSCAR technique is the sharp decrease of the MRFM relax- 
ation time when the distance between the ferromagnetic particle and the sample is below 
1 micrometer [1]. Thus, an important problem of the MRFM technique is to understand 
the nature of the spin relaxation in the submicron region and to find a way to reduce the 
relaxation rate. One of the possible mechanisms involved in MRFM relaxation in the submi- 
cron region is thermal currents in the metallic ferromagnetic particle [1]. Another possible 
mechanism is thermal excitation of the high frequency modes of the cantilever vibrations 



In this paper, we explore this second mechanism. We show that the thermal excitations 
of the high frequency modes may provide a noticeable contribution to the experimentally 
observed spin relaxation of the MRFM signal. We also show that the use of a nonuniform 
cantilever can reduce the spin relaxation rate caused by this mechanism. 



[5]. 
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FIG. 1. The geometry of the OSCAR technique. 
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II. MODEL AND EQUATIONS OF MOTION 



The geometry of the problem is shown in Fig. 1. This geometry is essentially the same as 
that considered in [2]. However, now we consider not just a single magnetic moment below 
the cantilever tip, but a uniform distribution of the magnetic moments across the sample 
which corresponds to the experimental conditions in [1]. The ferromagnetic particle attached 
to the cantilever tip effectively interacts with all magnetic moments inside the resonant slice 
where the Larmor frequency matches the frequency of the microwave magnetic field, B\. 
(The effective thickness of the resonant slice equals twice the amplitude of the cantilever 
vibrations.) 

The tip vibrations associated with the fundamental mode of the cantilever can be de- 
scribed as a motion of an oscillator with a quality factor Q: 



Here we used the dimensionless quantities: r = u c t, where oo c is the "unperturbed" cantilever 
frequency of the fundamental mode; z c is the z-coordinate of the center of ferromagnetic 
particle in the units of the cantilever amplitude z ; Xk, Uk, %k are the coordinates of the fc-th 
magnetic moment of the sample in the same units; Zk = Zk — z c ; r\ = x\ + y\ + z\. /i is 
the permeability of free space; jlk is the fc-th magnetic moment in units of its magnitude \x 
(which is the same for all magnetic moments); k c is the cantilever spring constant, m is the 
magnetic moment of the ferromagnetic particle. The function f{r) describes the action of 
the magnetic moments in the resonant slice on the cantilever. 

The motion of the /c-th macroscopic magnetic moment in the rotating reference frame is 
given by 



z c + z c + z c /Q = /(r), 



(1) 



f( T ) = YsVkHkz, 
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In Eqs. (2) the following notation was used: 



A fc = {jB ext -u)/u) c + 



fi ym 3g 

47T U c Zq 



(3) 



£ = -fB 1 /uj c . 

Here 7 is the electron gyromagnetic ratio in the sample; B ext is the external static magnetic 
field; B\ and u are the amplitude and frequency of the rotating rf magnetic field. Note, that 



The upper and the lower boundaries of the resonant slice are determined by the equation 

Afc = 0, at z c = ±1. 

In our numerical experiments, the magnetic moments have been distributed uniformly inside 
the resonant slice. We used the following initial conditions 



i.e. the magnetic moments are oriented approximately along the effective field in the rotating 
frame. To model the feedback technique in the OSCAR MRFM, our computer algorithm 
increased the value of z c to 1 every time the cantilever passed the upper point. The period 
of the cantilever oscillations was determined as the time interval between the instants of the 
z c maximum values. 

To model the thermal excitations of the high frequency cantilever modes we use the 
following approach. The cantilever energy E can be written as a sum of the energies of 
all vibrational modes. The high frequency modes can be described approximately by the 
formula [6] 



A k is the ^-component of the rotating-frame effective field B^^ (in units oo c /l), an d £ is the 
re-component of the effective field in the same units, = (e, 0, A k ). 
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[n(n - 0.5)] 2 , n > 1, 
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where t c is the thickness of the cantilever, / is its length, Y is Young's modulus, and p is the 
density of the cantilever. Combining the expressions 

= k B T, ujI = 4k c /m c , (6) 

where m c and k c are the mass and the spring constant of the cantilever, we estimate the 
thermal amplitude of the n-th mode to be 

1 / k T\ ^ 2 

an= VT\2^k~) '^« = W «M- (7) 

The amplitude of the cantilever tip oscillations, b n , for any high frequency mode, n, is twice 
the amplitude of the mode [6]: 

b n = 2a n = (2k B T/k c n 2 n ) 1 / 2 . 

To describe the influence of the noise on the spin dynamics we replace the coordinate z c 
with z c + 5z c in the expression for Zk and fk in (3), where 

fo c = £ — cos(ft n T + *n), (8) 

n z 

and ^ n being a random phase. 



III. RESULTS OF NUMERICAL SIMULATIONS 

We solved the system of equations (1), (2) changing z c — > z c + 5z c in Eqs. (2) according 
to (8). (We did not take into account the influence of the thermal noise on r/ k in Eqs. 
(1). As our simulations demonstrate, the influence of the thermal noise does not cause a 
significant direct contribution to the cantilever vibrations through the parameter r/k, but 
causes a dephasing of magnetic moments in the resonant slice.) The number of magnetic 
moments in the resonant slice was 200 (1 < k < 200). Our simulations show that the results 
do not change significantly when this number is increased to 400. 

The parameters in our numerical simulations were taken from experiment [1]: 

B ext = 140 mT, k c = 0.014 N/m, uj c /2tt = 21.4 kHz, z = 28 nm, (9) 
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Q = 2 x 10 4 , m = 1.5 x 1(T 12 J/T. 

The distance between the bottom of the ferromagnetic particle and the surface of the sample, 
di, (see Fig. 1) was d\ = 700 nm, and the radius of the ferromagnetic particle was equal to 
di. The average magnetization of the sample was 0.89 A/m. The rotating r/ field amplitude, 
Bi, in experiments was no less than 0.15 mT. The cantilever temperature, T, did not exceed 
20 K. (For these values of parameters the value of d<i in Fig. 1 was found to be di = 875 nm, 
and the value of the magnetic field gradient at the center of the resonant slice was 1.4 x 10 5 
T/m.) 

We have studied the decay of the OSCAR signal, AT/AT , where AT is the shift of 
the period of the cantilever vibrations due to the influence of the magnetic moments of the 
sample, and AT is the initial value of AT. In our model, the decay of the OSCAR signal 
is associated with cantilever noise near the Rabi frequency uor = ^B\ which is the smallest 
resonant frequency in the rotating reference frame. 

High frequency cantilever modes may generate significant noise near the Rabi frequency 
causing a deviation of the magnetic moments from the direction of the effective field and a 
nonuniform dephasing between magnetic moments. 

We write the amplitude of the cantilever tip thermal vibrations associated with the ra-th 
mode as b n = b R e/Q n , where 

b R = E- l (2k B T/k c ) 1 ' 2 . 

The phases *ff n in (8) were changed randomly between and 2n with random time 
intervals, r*, between two successive changes of the phase. In particular, we put r$ = Nto, 
where To takes random values between 2n/1.2e and 2n/0.8e, and A is a free parameter of 
the model. 

Our simulations show that the decay of the OSCAR signal is almost independent of N 
for N < 1000. (See Fig. 2.) We have found that for N < 1000 the signal decay can be 
approximately described by an exponential function with relaxation time, r m (Fig. 2). 
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FIG. 2. Decay of the OSCAR signal for different values of the parameter N. Curves 1-5 corre- 
spond to N = 10 5 , 10 3 , 100, 10 and 2. The amplitudes: bn = 5pm, z = 15nm, and e = 390. 

The relaxation time quickly decreases with an increase of the characteristic noise am- 
plitude at the Rabi frequency, b R which, in turn, is proportional to VT/e. Taking the 
minimum value of e (e — 195, i.e. B 1 = 0.15mT) and the maximum value of T (T = 20K) 
from the range of parameters [1], we obtain the maximum value of 6r, b R = 1pm. The 
corresponding minimum value of the relaxation time is approximately r m « 570ms. The 
experimental value found in [1], r m = 68ms, is smaller but of the same order of magnitude 
as our minimum theoretical value. Thus, we conclude that the thermal vibrations of the 
high frequency modes provide a noticeable contribution to the spin relaxation. 

To study the spin relaxation we also considered values of parameters which provide 
relatively small relaxation time. This allowed us to reduce the simulation time and to 
determine the characteristic scaling properties of the relaxation process. As an example, 
Fig. 3 shows the decay of the OSCAR signal for various values of the cantilever amplitudes. 
One can see that the relaxation time r m decreases when the cantilever amplitude decreases. 
(See curves a and c in Fig. 3.) Note that if we take into consideration spins near the resonant 
slice, the value of r m slightly increases. (Compare curves b and c in Fig. 3.) 
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FIG. 3. Decay of the OSCAR signal for various values of the cantilever amplitude zq; b^ = 1pm, 
and e = 390. Curve a corresponds to zq = 15nm; curves b and c correspond to zq = 7.5nm. For 
curve b we took into consideration spins close to the resonant slice. 
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FIG. 4. Decay of the OSCAR signal for zq = 28nm and various values of the temperature T and 
rotating magnetic field Bi. Curves 1-3 correspond to e = 390 (Bi = 0.3mT) and T = 20K, 40K 
and 80K. Curve 4 corresponds to e = 195 {B\ = 0.15mT) and T = 80K. The relaxation time in 
curves 1-4 is 1500ms, 700ms, 310ms, and 145ms. 
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Fig. 4 shows the decay of the OSCAR signal for various values of the temperature and 
rotating magnetic field. One can observe the expected decrease of the relaxation time r m 
with an increase of the temperature or decrease of the rotating magnetic field. 

Fig. 5 shows the decay of the OSCAR signal for various numbers of high frequency modes 
taken into consideration. One can see that the relaxation time r m slightly decreases when 
the number of high frequency cantilever modes increases from 2 to 22. A further increase of 
the number of modes does not essentially change the decay rate. 

Our simulations show that the OSCAR signal generally consists of two parts: (a) a 
regular part which we studied in our paper and (b) a random part which could be observed 
after the decay of the regular signal. Numerical analysis of the random part of the OSCAR 
signal requires large computational times, and will be studied elsewhere. 




FIG. 5. Decay of the OSCAR signal for 6r = 5pm, zq = 15nm, and various numbers of the 
high frequency modes taken into consideration. The lowest of high frequency cantilever modes is 
the mode with the frequency closest to the Rabi frequency 8.4MHz (e = 390). Curves a, b, and c 
correspond to 22, 3, and 2 high frequency cantilever modes including the lowest one. 
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IV. REDUCTION OF THE RELAXATION RATE 



In this section we discuss how to reduce the spin relaxation rate caused by thermal 
excitations of the high frequency cantilever modes. For this purpose we consider a cantilever 
with nonuniform thickness. 

The mechanical energy of the cantilever can be written as 



E — - 
2 



1 fi ' ^ 2 



o 



dx. (10) 



dt J \dx 2 / 

Here S is the cross-sectional area, / = wtl/12, and w is the width of the cantilever. The 
equation for the cantilever motion without an external force and damping is given by 

d 2 z d 2 ( d 2 z\ 



Using a new variable, £ = zyS(x), we obtain the equation of motion in the form 

d 2 C_ Y &_( ^d 2 ( \ n9 , 

P dt 2 ~ ^dx^dx 2 ^)- {L2) 

The general solution of the equation (12) can be written as 

C(x, t) = a n f n (x) exp(-iu n t). (13) 

n 

The eigenfunctions f n (x) are normalized in the following way: 

/ S(x)f n (x)f m (x)dx = V5 nm , (14) 
J o 

where V — Jq S{x)dx is the cantilever volume. The equation for the eigenfunctions is 

= -^=^( I ( X ^ J T^=)- ( 15 ) 

y/S(x) dX \ dx y/S(x)J 

We express the solution of Eq. (15) in the form 

fn(x) = c kn<fk(x), (16) 

fc=l 

where <Pk(%) is the fc-th eigenfunction for a cantilever with uniform cross-sectional area [6], 
satisfying the equation: 

SpulVn = YI^ (17) 
10 



/ 

Jo 



(Pn(x)<fm(x)dx = 15, 



nm- 



Taking k max = 70, we use a standard procedure to find the coefficients Ck n and frequencies 
uj n . Finally, using the same approach as for the uniform cross-sectional area we can find 
thermal vibrations of the cantilever tip, 5z c , in (8) with new values of Q n and b n . 

We found that it is possible to reduce the relaxation rate using a nonuniform thickness 
of the cantilever. (See also [5].) We considered a cantilever whose thickness increases near 
the location of the ferromagnetic particle: 



For the value a/l ~ 0.01, we obtained an almost tenfold increase in the relaxation time. 

We also considered the case when the ferromagnetic particle was located away from the 
tip of the cantilever. We found that in this case the increase of the cantilever thickness near 
the location of the ferromagnetic particle causes a significant increase of the relaxation time. 
The details of these results will be published elsewhere. 



We simulated the spin relaxation in OSCAR MRFM technique caused by the thermal 
excitation of the high frequency cantilever modes. Using a range of experimental parameters 
[1] we obtained the minimal relaxation time 570 ms which is greater than but of the same 
order of magnitude as the experimental value 68 ms. We demonstrated that an increase 
of the cantilever thickness at the location of the ferromagnetic particle can significantly 
increase the relaxation time. 

Finally, we would like to mention two assumptions of our computational model: We 
assumed 1) that the mass of the ferromagnetic particle was small compared to the cantilever 
mass, and 2) that the ferromagnetic particle has a spherical shape. We are now simulating 
a modified cantilever-spin model which is not restricted by these two assumptions. 




(18) 



V. CONCLUSION 
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